Nanowire formation on sputter eroded surfaces 
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Rotated ripple structures (RRS) on sputter eroded surfaces are potential candidates for nanoscale 
wire fabrication. We show that the necessary condition for RRS formation is that the width of the 
collision cascade in the longitudinal direction has to be larger than that in the transverse direction, 
which can be achieved by using high energy ion beams. By calculating the structure factor for the 
RRS we find that they are more regular and their amplitude is more enhanced compared to the 
much studied ripple structure forming in the linear regime of sputter erosion. 
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The fabrication of nanoscale surface structures such as 
quantum dots (QDs) and quantum wires (QWs), have at- 
tracted considerable attention due to their applications 
in optical and electronic devices These nanostruc- 
tures form thanks to various self-assembled mechanisms, 
induced by the combined effect of strain and growth ki- 
netics. Yet, the strained nanostructures obtained by 
these methods have a size distribution wider than re- 
quired by applications, and display random alignment. 
Lithographic methods |^ are often considered prime can- 
didates to overcome these shortcomings, but their limited 
resolution offers further challenges. Consequently, there 
is continued high demand for alternative methods that 
would allow low cost and efficient mass fabrication of 
nanoscale surface structures. In the light of these tech- 
nological and scientific driving forces, the recent demon- 
stration by Facsko et al. that low-energy (40 eV ~ 1.8 
keV) normal incident Ar+ sputtering on GaSb (100) sur- 
faces leads to nanoscale islands which display remarkably 
good hexagonal ordering and have a uniform size distri- 
bution, has captured the interest of the scientific com- 
munity |,|. 

It is known that the morphological evolution of a sput- 
ter eroded surface is well approximated by the noisy non- 
linear Kuramoto-Sivashinsky (KS) equation. 



dth = v^dlh + Vydyh - D^^d^h - Dyydyh - D^ydldyh 

<i^,y,t), (1) 



where and Vy are the effective surface tensions gen- 
erated by the erosion process; D^^, Dyy, and D^y are 
the ion induced effective diffusion constants; and Aj, 
describe the tilt-dependent erosion rates in each direc- 
tion; and ^{x,y,t) is an uncorrelated white noise with 
zero mean, mimicking the randomness resulting from the 
stochastic nature of ion arrival to the surface |^,^. At 
low temperatures all the coefficients in Eq. (|l|) depend on 
experimental parameters such as the ion beam flux /, the 
ion beam energy e, and the incidence angle of ion beam 
6 [0, while at high temperatures D depends on surface 
temperature. For the sputter erosion process v < and 
D > 0, while the signs of A^ and Ay vary depending on 
the incident angle of the ion beam. 

Recently numerical simulations have shown that there 



is a clear separation of the linear and nonlinear regimes 
in time Up to a crossover time ti the surface is 
eroded as if the nonlinear terms would be completely 
absent, following the predictions of the linear theory |^ 
(i.e. A3; = Ay = in Eq. (|l|)). After ri, however, the 
nonlinear terms with coefficients A^; and Ay take over and 
completely determine the surface morphology. The tran- 
sition from the linear to the nonlinear regime can be seen 
either by monitoring the surface width (which is propor- 
tional to the ripple amplitude) or the erosion velocity. 

In the nonlinear regime the case Aa;Ay < is in par- 
ticular interesting. The surface morphology in this case 
exhibits another transition from kinetic roughening to 
a rotated ripple structures (RRS) at a second crossover 
time T2 {t2 > Ti), as first predicted by Rost and Krug 
JlOt , and observed numerically in Moreover, it was 
found that the RRS is straighter than the ripple pattern 
forming in the linear regime, making it a potential candi- 
date for nanowire fabrication. In this Letter we investi- 
gate the necessary conditions for the formation of RRSs 
and the impact of the various experimentally controllable 
parameters on the morphology of the RRS. We find that 
at low temperatures the RRSs form when the longitudi- 
nal width (cr) of the damage cascade generated by the ion 
beam is larger than the transverse width (/i), which can 
be achieved for high energy (e). When the energy is too 
high, however, the crossover time T2 becomes too long. 
Thus, we predict that the RRS structure could be ob- 
tained in a moderate range of the incident energy, using 
an appropriate exposure time and low temperatures. 

Rotated ripple structures form when A^jAy < and 
the rotation angle is given by (j>c = tan""'^ ^J—Xx/Xy- As 
shown in Fig. 1, the angle 0c increases with the ratio 
= a/ 11, but decreases with the incident angle 9, where 
a denotes the penetration depth of the ion beam and 
depends on the ion beam energy e. The KS equation in 
the rotated frame can be written in the same form as 
Eq. (P except that the coefficients D, and A are re- 
placed by v' , D', and A', which are functions of those in 
the original frame and the angle (jjc- In the rotated frame 
one of the coefficients of the nonlinear terms, say X'^, , is 
equal to zero, and the other is given by X'^, = X^ + Xy, 
where {x',y') represents the coordinates in the rotated 
frame. Since A^/ vanishes, the dynamic equation in the 
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x' direction becomes linear. Consequently, the ripple pat- 
tern is along the x' direction as long as (i) v'^, < 0; (ii) 
ZJ^/j./ > 0; and (iii) Xx\ < 0. Therefore the conditions 
(i), (ii), and (iii) are the necessary conditions for the for- 
mation of the RRS. 

We investigate the satisfiability of these conditions in 
the parameter space {0,a^ = a/fi) for different values of 
Oct = a/ a. We find that when > 1 and = 1 or 
> 2 and flo- = 2, the RRS can form in the region de- 
picted in Fig. 2. For < 1 given — 1, the shaded 
region satisfying (i) - (iii) scarcely exists, so that the for- 
mation of RRS is less likely. That means the RRSs are 
expected to form when the longitudinal width a is larger 
than the transverse width fi, that is, cr > /i. Recent 
experimental results indicate that for graphite surfaces 
a depends on e, while fi is independent of e for large e 
(2 ~ 50 keV) [|ll|. Therefore, the a > fi condition can 
be met when the energy of the incident ion beam is high 
enough. Thus in order to obtain the RRS experimentally 
it is desirable to use high energy ion beam with an ap- 
propriate choice of the incident angle (Fig. 2). However, 
the use of a high energy ion beam increases T2 rapidly, 
requiring a longer exposure time. 

Since the nonlinear term disappears in the x' direction 
the surface in this direction is driven by a linear insta- 
bility. The amplitude of the RRS grows exponentially 
with time until the nonlinear term in the y' direction 
becomes effective, after which the amplitude saturates. 
Meanwhile, the surface in the y' direction displays ki- 
netic roughening due to the presence of the nonlinear 
term Ay (9y//i)^/2, so that the roughness in the y' direc- 
tion is considerably reduced compared with the rough- 
ness in the x' direction. Therefore the RRS develops a 
rough morphology in the x' direction, while it is rela- 
tively smooth in the y' direction, the end configuration 
resembling a V-shaped wire pattern, as shown in Fig. 3 
[ p2t . The RRS formed in the nonlinear regime is compa- 
rable with the ripple pattern formed in the linear regime, 
where there are modulations in both directions, and the 
roughness in each direction is almost of the same order. 

We also examined the structure factor 

^^"^^^ / (|^exp(zq.r)i/(r), (2) 
where H{r) is the height-height correlation function, 

^('■'*) = ^EMr' + r,t)Mr',<)^ - (^E ' 

(3) 

averaged over different configurations We find that 
the structure factor exhibits a peak at (qx,c, Qy,c), corre- 
sponding to (g^/ 0) in the rotated coordinaates (Fig. 4). 
The fact that q'^, ^ = confirms that the RRS is straight 
in the y' direction. Moreover, the amplitude of the struc- 
ture factor for the RRS is much larger compared to the 



ripple formed in the linear regime, as shown in the inset 
of Fig. 4. Accordingly, the RRS could be a good candi- 
date for the fabrication of nanowires. 

Applying the linear instability theory in the x' di- 
rection we obtain the wavelength of the RRSs £'^, — 
2t^^2D[j.i,j,, l\v'^,\. The wavelength depends on the pen- 
etration depth linearly, i.e. i'^i ^ a, independent of the 
ion flux /. In general, the penetration depth depends on 
the incident ion energy e as a ~ e^™ , where m « 1/4 was 
obtained recently for the nanoscalc dot structure in low 
energy sputtering fl^ . We also examined in function 
of the ratio and incident angle 6, observing a mono- 
tonically increasing and decreasing behavior, respectively 
(Fig. 5). 

In the high temperature limit the wavelength of the 
RRS depends on the ion energy as ~ e~^^^, and 
on the ion flux as ~ /~^^^, and on the temperature as 
~ exp(— e/4A;T), a dependence similar to that observed 
during the formation of nanoscale dots ||l^ . 

In summary, we examined the necessary conditions 
for the formation of RRSs, potential candidates for 
nanowires for electron transport. We predict that the 
RRS can be generated under high energy ion beam, in 
contrast with the formation of nanoscale dots structure 
occurring during low energy ion beam sputtering. Since 
high energy ion beam causes longer crossover time for the 
formation of the RRS, exposure time of the ion beam has 
to be adjusted to obtain the desired RRS. 
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FIG. 3. Surface morphology of the RRS, as generated by 
numerical simulations, with ay, = 1.3, = 1 and 6 = 43.56 
degree. 




FIG. 1. The rotation angle < 



as a function of a,j at i 



50 

degree. The inset shows 0c versus 9 for = 1.5. The rota- 
tion angle increases as the ratio increases, but decreases 
with the incident angle 6. 
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FIG. 2. The shaded region in the parameter space (6',a^) 
for tto- = 1 (inset: for Oa = 2) corresponds to the region where 
the RRS can form. 
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FIG. 4. The amplitude of the structure factor jS'(q')| for 
the RRS shown in Fig. 3. The peak of the structure factor is at 
'^l'x',c — 0-33 and q'y, ,, = 0, implying that the wire structure 
is straight along the y' axis. The inset shows the comparison 
between the amplitudes of the structure factor for the RRS 
(solid line) and for the linear ripple (dotted line) , implying the 
amplitude of the RRS is about a factor of 10^ larger compared 
to that of the ripple structure formed in the linear regime. 




FIG. 5. The wavelength of tl^e RRS in units of a as a, 
showing function of (6 = 50 degree). Inset: The wave- 
^ length versus ^ at a^^ = 1-5. 



